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Abstract. We use the possible Cepheid-like luminosity estimator for the long-duration gamma- ray bursts (GRBs) developed 
by Reichart et al. (2000) to estimate the intrinsic luminosity, and thus the redshift, of 907 long-duration GRBs from the 
BATSE 4B catalog. We describe a method based on Bayesian inference which allows us to infer the intrinsic GRB burst rate 
as a function of redshift for bursts with estimated intrinsic luminosities and redshifts. We apply this method to the above 
sample of long-duration GRBs, and present some preliminary results. 



INTRODUCTION 

There is increasing evidence that gamma-ray bursts 
(GRBs) are due to the collapse of massive stars (see, e.g., 
[1] for a discussion of this evidence). If GRBs are indeed 
related to the collapse of massive stars, one expects the 
GRB rate to be roughly proportional to the star-formation 
rate (SFR). However, the observed redshift distribution 
of GRBs differs noticeably from that of the SFR: the ob- 
served GRB redshift distribution peaks at z w 1 and few 
bursts are observed beyond z 1 .5, while the SFR peaks 
at z « 2 and 10-40% of stars are thought to form beyond 
z = 5 (see, e.g., [2,3]). 

However, observational selection effects play a key 
role in determining the observed redshift distribution of 
GRBs. Among these selection effects are the efficiencies 
with which burst redshifts can be determined by spectro- 
scopic observations of the burst afterglow and/or the host 
galaxy of the burst. In addition, both of these methods of 
determining redshift require the identification of an op- 
tical afterglow of the GRB, and the detectability of the 
optical afterglow may be a function of redshift. 

In this paper, we take a different approach. We use 
the possible Cepheid-like luminosity estimator for the 
long-duration GRBs developed by Reichart et al. [4] 
to estimate the intrinsic luminosity, and thus the red- 
shift, of 907 long-duration GRBs from the BATSE 4B 
catalog [5]. This approach is free of the important and 
difficult-to-quantify observational selection effects that 
affect the redshift distribution of the GRBs with known 
redshifts. We assume a very general model for the GRB 
rate and a truncated power-law model for the intrin- 
sic GRB isotropic-equivalent photon luminosity distri- 
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FIGURE 1. The variabilities V and isotropic-equivalent peak 
photon energy luminosities L of the 19 bursts for which some 
redshift information exists. The solid and dotted lines mare the 
center and l-CT widths of the best-fit model of these bursts in 
the (log L, log V)-plane. 



bution (in this work we assume that the amplitude and 
the power-law index of the isotropic-equivalent photon 
luminosity distribution does not evolve; we will relax 
this assumption in future work). Adopting a Bayesian ap- 
proach, we calculate the likelihood of the data given the 
model, and convert it to a posterior distribution on the 
model parameters. 
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FIGURE 2. Left panel: distribution of the 907 long-duration BATSE bursts in the (log z, log L)-plane, as determined by our 
luminosity estimator. The solid and dashed curves represent the 90% and 10% detection thresholds of BATSE. The horizontal 
dotted lines show the values of and Lniax adopted in this study. The diamonds show the locations of the 19 GRBs with known 
redshifts used to calibrate the luminosity estimator. Right panel: a Monte Carlo simulation of the expected distribution of f« 900 
bursts with peak photon fluxes P above the BATSE 90% detection threshold, assuming the maximum likeUhood best-fit parameters 
for the model (see text). 



LUMINOSITY ESTIMATOR 

A possible Cepheid-like luminosity estimator for 
gamma-ray bursts has been suggested by Ramirez-Ruiz 
and Fenimore [6] and developed further by Reichart et 
al. [4,7,8]. These authors have shown that there exists a 
correlation between a measure V of the variability of the 
burst time history and the intrinsic isotropic-equivalent 
peak photon energy luminosity L of the burst for the 
19 GRBs for which some redshift information (either 
a redshift measurement or a redshift hmit) exists. We 
show this correlation in Figure 1. 

We apply this luminosity estimator to 907 bursts from 
the BATSE 4B catalog [5] that have durations Tg^ > 10 
sec. We show the resulting redshift-luminosity distribu- 
tion for these bursts in Figure 2 (see also [7] and [8]). 



MODEL 

In this work, we assume that the rate of GRBs per unit 
redshift and luminosity is given by a separable func- 
tion of redshift z and intrinsic isotropic-equivalent pho- 
ton number luminosity L. That is, we assume that the 
intrinsic photon number luminosity distribution of GRBs 
does not evolve with redshift; we will relax this assump- 
tion in future work. Then the rate of GRBs per unit red- 



shift and luminosity can be written as [9] 
dN 



dzdL^ 
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where 

P(z) =^GRB(2;P.e) X {l+z)-'x4Kr{zf{dr/dz) (2) 

is the rate of GRBs that occur at redshift z, and r{z) is the 
comoving distance to the source. 

We adopt a phenomenological model for the rate of 
GRBs that occur at redshift z per unit comoving volume 
of the form proposed by Rowan-Robinson for the star 
formation rate [10]. In this model, the GRB rate is given 
by 
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where P, Q and are model parameters, and t{z) is the 
time since the Big Bang corresponding to the redshift z- 
We take the intrinsic photon luminosity distribution of 
GRBs to be a truncated power-law. 
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Thus the model has six parameters: P, Q, Rq, j3, L^^^, 
andLmax- 
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FIGURE 3. Left panel: the best-fit point in the {P, 2)-plane at the maximum likelihood best-fit values of the normalization 
Rq = 3.46 of the GRB rate and the power-law index ^ = 1 .52 of the GRB luminosity distribution. The contours shown correspond 
to Alog^ = 150. Right panel: a comparison of the expected and observed number of bursts with P above the BATSE 90% detection 
threshold for the maximum likelihood best-fit model. 



In this work we assume a flat universe with Q.^ = 0.3 
and Q.^ = 0.7. Furthermore, we fix L^^, and Lmax to 
be constants (see Figure 2), reducing the number of free 
parameters in our model to four. 



function, to be 8{P — P^^^), where =0.4 photons cm 
sec"i is the 90% BATSE detection threshold. The prod- 
uct of the event likelihoods jx{zi,Pi) runs over the 907 
BATSE bursts considered in this study. 



LIKELIHOOD FUNCTION 

The likelihood function for GRBs that are Poisson in 
time is given by [9] 

^ = exp|-| dzdPii{z,P)e{z,P)^fY{Zi,Pi) , (5) 
where the event likelihood, 

= p(z,)x/(47rr(z,)2(l+z,)«/^-) 

x47trizfil+Zi)", (6) 

is the expected number of events observed within dzdP 
of redshift z, and peak photon number flux P^. The quan- 
tity a is the burst spectral index, which we take to be 
equal to 1.0, a value that is typical of GRBs [11]. By an 
application of Bayes' Theorem, we regard ^ as an (un- 
normalized) probability distribution on the model param- 
eters. We estimate the best-fit parameters in our model by 
maximizing this likelihood function over the parameter 
space. 

We take the BATSE observing efficiency e{z,P), 
which appears in the exponential factor in the likeUhood 



RESULTS 

Each (z, L)-point in the left panel of Figure 2 corresponds 
to the maximum value of a probability distribution for 
that particular burst. Thus each point has statistical and 
systematic errors that are associated with it and that are 
not displayed in the figure. Since the estimate of the red- 
shift z is calculated from the estimate of the peak photon 
luminosity L (and the peak observed photon number flux 
P), the errors in z and L are completely correlated; that 
is, the uncertainty lies along curves of constant P (i.e., 
they are diagonals) in the {z, L)-plane. Furthermore, these 
uncertainties are not symmetric; they are skewed toward 
low L at low z and toward high L at high z [7,8] . 

In this preliminary work, we neglect these statistical 
and systematic errors, and use only the best-fit values in 
our analysis. Calculating the likelihood of this data, given 
the model, and converting it to a posterior distribution on 
the model parameters, we find the maximum fikefihood 
parameters for this data set to be: P = —0.70, Q = —0.55, 
Rq = 3.46 and j3 = 1.52. The right panel of Figure 2 
shows a Monte Carlo simulation of the distribution of 
w 900 bursts whose peak photon fluxes P are above 
the BATSE 90% detection threshold for the maximum 
likelihood best-fit parameters. 
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FIGURE 4. Left panel: the best-fit GRB rate per unit comoving volume (solid line). Right panel: the GRB rate per unit redshift z. 
Estimates of the star formation rate as a function of redshift z made by Madau et al. (1998) and Rowan-Robinson (1999) are shown 
in both panels for comparison, as is the no evolution with redshift z model {P = Q = 0). We emphasize that we have not taken into 
account the statistical and systematic errors in the redshifts z and intrinsic peak photon luminosities L derived from the variabiUty 
measure V, and therefore cannot quote meaningful confidence regions for our best-fit parameters. 



The left panel of Figure 3 shows the best-fit point 
and contours of Alog^ = 150 in the (P,2)-plane at 
the maximum Hkelihood best-fit values of the normal- 
ization Rq = 3.46 of the GRB rate and the power- law 
index J3 = 1 .52 of the GRB luminosity distribution. The 
left panel of Figure 3 thus represents a two-dimensional 
slice through the four-dimensional parameter space of 
the model. The contours shown in Figure 3 do not corre- 
spond to credible regions because we have not taken into 
account the statistical and systematic errors in the red- 
shifts z and intrinsic peak photon luminosities L derived 
from the variabihty measure V. However, taking the pre- 
liminary results at face value, the best-fit model appears 
to be consistent with no evolution but to be inconsistent 
with the Rowan-Robinson [10] and the Madau et al. [12] 
expressions for the star formation rate as a function of z- 

The right panel of Figure 3 shows a comparison of the 
observed and expected numbers of bursts with P above 
the BATSE 90% detection threshold for the maximum 
likelihood best-fit model. These distributions thus corre- 
spond to projections onto the z-axis of the burst distribu- 
tions in the left and right panels of Figure 2, respectively. 

The left panel of Figure 4 shows the best-fit GRB rate 
per unit comoving volume, while the right panel shows 
the GRB rate per unit z- Estimates of the star formation 
rate as a function of z made by Madau et al. (1998) and 
Rowan-Robinson (1999) are shown for comparison, as is 
the no evolution model (P = Q = Q). 

Applying the intrinsic peak luminosity estimator of 
Ramirez-Ruiz and Fenimore [5] to the brightest 220 
long-duration bursts in the BATSE 4B catalog [11], 



Lloyd-Ronning et al. [13] find that the GRB rate in- 
creases with increasing redshift (Rq^q °^ (1 +z)^ until 
z ^ 2 and (1 from z ^ 2 until z ^ 10), and that 
the mean luminosity of the GRBs also increases with 
increasing redshift (< L (1 + z)' "*^" -^). The results 
found in the present preliminary work differ somewhat 
from these conclusions. But we again emphasize that (1) 
the model that we have used in the present preliminary 
study does not allow for the evolution of L with z, and 
(2) we have not taken into account the statistical and sys- 
tematic errors in the redshifts z and intrinsic peak pho- 
ton luminosities L derived from the variability measure 
V, and therefore cannot quote meaningful confidence re- 
gions for our best-fit parameters. In future work, we will 
relax the assumption of no evolution and will take into 
account the statistical and systematic errors. 
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